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摘要 :对 径流 演变 规律 的 认识 和 研究 是 水 资源 得 到 有 效 利 用 的 前 提 。 密 云 水 库 是 北京 市 唯一 的 地 表 水 源 地 ;水 库 的 正常 运行 对 
北京 市 供水 安全 有 重要 的 意义 。2015 年 9 月 ,南水北调 的 水 开始 进入 密云 水 库 , 对 密云 水 库 上 游 径流 量 的 研究 可 以 使 南 水 北 
调 的 水 与 其 相互 配合 ,有 效 补充 密云 水 库 的 水 资源 ,研究 上 游 河流 的 径流 变化 趋势 很 有 必要 以 密云 水 库 白河 ( 张 家 坟 站 ) 和 
潮 河 (下 会 站 ) 水 文 站 1960 一 2014 年 的 月 径流 数据 为 基础 ,采用 年 内 径流 不 均匀 系数 、 滑 动 平均 法 、Mann-Kendall 突变 检验 法 、 
小 波 分 析 方 法 . 双 累 积 径流 曲线 法 分 析 了 密云 水 库 入 库 径 流 的 年 内 分 配 特 征 和 年 际 变化 特征 ;径流 变化 周期 趋势 及 降雨 和 人 类 
活动 对 径流 变化 的 影响 。 结 果 表 明 :密云 水 库 入 流 径 流量 主要 分 布 在 7 一 9 月 之 间 ; 册 全 年 径流 量 的 70% ;白河 流域 年 内 径流 变 
化 较为 剧烈 。50 年 来 白河 与 潮 河 径流 量 有 减少 趋势 ,递减 率 分 别 为 3.24m s™(10a)” 和 二.60m s""(10a)”。 白 河流 域 和 潮 河 
流域 径流 量 存在 周期 为 7a 的 变化 特征 。 利 用 累积 径流 曲线 分 析 降 雨 和 人 类 活动 对 密云 水 库 上 游 径 流量 的 影响 ,结果 表明 人 类 
活动 是 影响 密云 水 库 径 流 变 化 的 主要 原因 ,为 71.36% 一 91.39% ,其 中 白河 流域 受 人 类 活动 干扰 更 为 剧烈 。 
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Abstract: Studying thetrends and laws of runoff is fundamental to the proper use of water resources. The Miyun Reservoir is 
the only surface watér source in Beijing and is important to Beijing's water supply. In September 2015, water from South- 
North Water Diversion Project (SNWDP) was added to the reservoir. Therefore, studying the runoff of upstream in Miyun 
Reservoir ls -more necessary for the proper management of the reservoir than ever before. As the amount of the water from 
SNWDP shiould supplement the runoff of the upstream, to know the trend of the upstream 1s necessary. Accordingly, the 
study-analyzed the multi-temporal scale changes and abrupt change of the runoff, the influence factors of the runoff of Zhang 
Jiafen and Xiahui stations from 1960 to 2014. The methods of trend analysis, the moving average method, the Mann-Kendall 
meéthod, wavelet analysis, the flow duration curve, and the double mass curve were used. The results indicated that the 
runoff“in Miyun Basin mainly concentrated in July to September, accounting for 70% of the annual runoff, and runoff in the 
Baihe Basin changed more significantly than that of the Chaohe Basin. From 1960 to 2014, the amount of runoff slowly 
decreased at the rate of 3.24 and 1.60m’ s (10a)™ in Baihe and Chahe respectively, and the runoff in Baihe and Chaohe 
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Basin both showed change period of 7a. The results of double mass analysis indicated that the reduction of runoff in Miyun 
Basin can mainly be attributed to human activities (71.36 一 91.39% ) ，and the Baihe Basin was disturbed much greater by 


human activities. 
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河川 径流 的 变化 是 气候 地貌 .土壤 及 人 类 活动 共同 作用 的 结果 '"。 分 析 河 川 径 流量 的 时 间 变 化 特征 、 
周期 规律 以 及 影响 因素 ,对 了 解 水 资源 特性 ,从 而 更 好 的 开发 利用 水 资源 有 重要 的 意义 '"。 尤 其 我 国 北方 地 
区 ,水 资源 已 是 制约 社会 经 济 发 展 的 关键 因素 … 。 密 云 水 库 是 北京 市 目前 唯一 的 地 表 饮 用 水 源 地 ,上 游人 库 
径流 的 变化 对 其 水 量 有 重要 的 影响 ,分 析 评 价 密云 水 库 和 库 径 流 的 年 际 变化 及 周期 性 演变 趋势 识别 影响 径 
流 变化 的 因素 对 优化 密云 水 库 水 资源 调度 ,可 为 北京 市 饮用 水 保障 提供 有 利 的 科学 依据 。 

日 前 对 密云 水 库 上 游 流 域 展 开 的 研究 采用 的 方法 主要 有 统计 经 验 公式 法 “水文 模 型 法 ;分 为 集 总 式 模 
型 “和 分 布 式 模型 中 。Li 等 “采用 M-K 检验 的 方法 分 析 了 潮 河 流域 降雨 和 径流 的 变化 趋势 ,表明 人 类 活动 
对 径流 变化 产生 重要 的 影响 ; 孙 宁 等 利用 IHACRES 模型 分 析 了 土地 利用 变化 和 履 被 对 潮 河 流域 径流 的 影 
响 ,结果 表明 ,土地 利用 和 覆 被 变化 减少 了 流域 年 径流 ; 郭 军 庭 等 “利用 SWAT 模 型 分 析 了 济 河 流域 径流 对 
气候 变化 和 土地 利用 变化 的 啊 应 特征 ,结果 表明 ,气候 变化 对 年 径流 量 的 影响 较 大 。 

然而 这 些 研 究 将 研究 区 域 集中 在 潮 河流 域 ””, 对 白河 流域 径流 特征 的 研究 较 少 。 特 别 是 对 影响 密云 水 
库 径 流量 的 影响 因素 仍然 存在 争议 , Wang 等 ”认为 潮 白 河流 域 1980 年 前 后 径流 减少 的 主要 原因 是 人 类 活动 
(68% 一 70% ) ,而 气候 变化 是 次 要 原因 (30% 一 32%)。 因 此 ;人 类 活动 对 潮 河 流域 径流 减少 起 决定 因素 。 而 
Ma 等 ”认为 潮 白 河流 域 径流 下 降 的 原因 主要 是 气候 变化 (55% 二 51%) ,而 人 类 活动 对 径流 下 降 的 影响 只 有 
18%。 特 别 是 密云 水 库 担任 南水北调 进 京 水 量 调 蓄 的 重要 任务 ,通过 对 上 游 来 水 量 的 径流 特征 深入 了 解 , 对 
南水北调 入 库 流 量 的 科学 调度 有 重要 的 意义 。 因 此 。 有 必要 对 密云 水 库 白河 和 潮 河 流域 的 径流 变化 趋势 和 
影响 因素 进行 分 析 和 比较 ,这 对 高 效 和 合理 管理 密云 水 库 有 重要 作用 。 


1 研究 区 域 概况 


密云 水 库 (40°29.0' 一 40°30YY 116$50.0' 一 117°3.5'E) 位 于 北京 北部 ,建成 于 1960 年 ,是 以 综合 防洪 、 
供水 为 目的 的 大 型 水 利 五 程 ,总 库容 为 43.7 亿 m”"”。 白 河和 潮 河 是 密云 水 库 主 要 的 入 库 河流 (图 1) ,其 
中 ,白马 关 河 、 蛇 渔 川 河 是 白河 的 主要 支流 ; 红 门 川 河 、 清 水 河 、 安 达 木 河 、 凡 牛 河 汤 河 是 潮 河 的 主要 支流 。 安 
达 木 河上 建 有 遥 桥 峪 水 库 , 虹 和 牛 河中 游 建 有 半 城 子 水 库 。 流 域 控 制 流域 面积 13788km ,属于 半 干 旱 半 湿润 的 
大 陆 性 季风 气候 ;降水 主要 分 布 在 汛期 (7 一 9 月 ) , 占 总 降水 量 的 80%。 流域 的 主要 土壤 类 型 为 褐 土 和 棕 壤 ， 
主要 植被 覆盖 类 型 是 针 阔 叶 混 交 林 |。 


2 数据 来 产 和 研究 方法 


2.1 ， 数据 来 源 

径流 量 资 料 来 源 于 密云 水 库 管理 处 ,包含 白河 与 漳 河 的 流域 出 口 张 家 坟 和 下 会 两 个 水 文 站 的 月 径流 资 
料 ,径流 系列 长 度 为 1960 一 2014 年 。 
22 研究 方法 

本 文采 用 的 研究 方法 包括 不 均匀 系数 法 分 析 年 内 径流 量 的 变化 趋势 ,利用 滑动 平均 法 分 析 径 流量 的 年 际 
变化 趋势 ,采用 M-K 突变 分 析 法 确定 径流 突变 点 ,利用 Molert 小 波 分 析 50 年 来 密云 水 库 径 流 的 周期 性 变化 ， 
并 通过 对 周期 性 的 显著 性 检验 ,确定 主 周期 。 利 用 降水 -径流 双 累 积 曲 线 分 析 和 揭示 降水 和 人 类 活动 对 径流 
的 影响 阶段 性 变化 ,从 而 得 到 降水 变化 和 人 类 活动 变化 对 径流 的 贡献 ,具体 流程 见 图 2。 
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| miles 
po—< 图 例 
人 一 一 水 系 | |] 密云 水 库 流域 加 密云 水 库 
«TT 
CO 图 1 密云 水 库 位 置 图 
Fig.1 The location of Miyun Reservoir 


CI ”2.2.1 年 内 变化 
ES 径流 量 年 内 不 均匀 系数 可 以 通过 公式 (1) 一 公式 
S 志 (3) 进 行 计算 。 


I KS OO 
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= | 本 径流 基本 特征 
各 ] 七 
O 5= /> ER: -As (2) 径流 年 际 变化 趋势 
i=1 
1 12 
ey (3) MK 突变 分 析 
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式 中 ,Re 为 月 平均 径流 量 ( 亿 mm ) ,RR 为 第 i 月 的 径 
流量 ( 亿 An ) ,0 为 均 方差 系数 
2.2.2/ 年 际 变化 

水 文 序列 的 随机 波动 性 使 得 难以 从 序列 本 身 判 断 
出 其 变化 趋势 。 滑 动 平均 法 可 以 消除 波动 的 影响 ,通过 
对 原 序列 的 光滑 化 处 理 , 即 对 序列 的 几 个 前 期 值 和 后 期 
值 取 平 均 ,得 到 新 的 序列 ,通过 新 序列 判断 水 文 变化 的 


径流 周期 性 
分 析 


径流 影响 因素 和 
> 双 累积 径流 曲线 


趋势 性 具有 直观 和 明显 的 优势 。 利 用 公式 (4) 进行 图 2 密云 水 库 径流 特征 分 析 流程 图 
计算 Fig.2 The flow chart of the analysis method of Miyun runoff 
] < 
vi (4) 


式 中 ,为 滑动 长 度 ,j=(1,2,3,…,n-k+1) 
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2.2.3 ”突变 分 析 -Mann-Kendall 非 参 数 统计 检验 

Mann-Kendall 非 参 数 统计 检验 方法 ,可 进行 径流 量 趋势 的 检验 ,其 优点 在 于 样本 不 必 遵 循 某 一 特定 规律 ， 
也 不 受 个 别 异常 值 干扰 ,能 客观 表征 样本 序列 的 整体 变化 趋势 " 。 

设 某 一 时 间 序 列 X 具 有 nn 个 样本 量 ,构造 该 时 间 序 列 的 秩序 列 . 


h 
5 2 (KS 1 2 (5) 
i=1 
3 
I 
六 下 人 三 人 2 (6) 
0 wx; < 和 


秩序 列 5, 表 示 第 i 时 刻 序列 值 大 于 j 时 刻 序列 值 的 个 数 之 和 。 当 x ,x ,…,*, 相 互 独立 且 连 续 同 分 布 时 ， 
定义 统计 量 ， 


[90 BCD) 
DR C7) 
a Vor 
2 (8) 
Re 2 (9) 


式 中 ,k=2,…,n; n 为 样本 数 ;UF 为 M-K 统 计量 , 且 VE = 0 且 WH, 服从 正 态 分 布 ;E(5,) 和 Var(S;) 分 别 是 
累计 数 5, 的 均值 和 方差 。 

取 显 若水 平 a, 奉 | UF,| > U。( 其 中 ,Uoos=+1.96,Woo<+2:57D95 则 说 明 序列 趋势 线 显著 。UF, 为 正 时 , 表 
明 序 列 呈 上 升 趋势 ,否则 呈 下 降 趋 势 , 且 其 绝对 值 越 大 表明 序列 变化 趋势 越 显 著 。 将 时 间 序 列 逆序 ,重复 M-K 
的 计算 过 程 ,得 到 逆序 列 M-K 统计 值 UB,， 当 统计 量 序列 曲线 出 现 交 点 , 且 交 点 在 临界 线 之 间 时 ,此 时 交点 所 
在 的 年 份 即 为 可 能 发 生 突变 的 时 间 ” 。 
2.2.4 周期 性 变化 

由 于 水 文 演变 时 序 是 多 时 间 尺 度 的 持续 过 程 ,而 传统 的 分 析 水 文 要 素 周 期 变化 特征 的 方法 对 揭示 其 
多 时 间 尺 度 特 征 存在 局 限 性 ,不 能 揭示 水 文 要 素 多 时 间 尺 度 的 特征 ,而 小 波 分 析 在 研究 不 同时 间 尺 度 变 换 
及 对 信号 的 目 适应 方面 有 具有 优势 因此 采用 小 波 分 析 的 方式 确定 密云 水 库 日 河 与 漳 河 入 库 径 流量 的 周期 性 
变化 。 主 要 的 小 波 函 数 有 "Daubechies 小 波 、Meyer 小 波 、Mexicanhat 小 波 、Symlet 小 波 、Morlet 小 波 等 :其 中 
Morlet 小 波 在 水 文 周期 性 序列 上 较为 常用 “ 。Morlet 小 波 函 数 为 复数 小 波 ,表达 式 如 下 "”: 

We (10) 

式 中 :e 为 常数 为 虚 部 :小 波 变换 是 对 于 基本 小 波 函数 做 位 移 后 ,再 对 不 同 尺度 时 域 下 的 需要 进行 分 析 的 
言 号 做 内 积 。 任 意 函 数 /(1) eL(R),L(R) 表 示 平 方 可 积 的 实数 空间 。f(1) 连续 小 波 变换 的 表达 式 为 公式 
(11) 


wa 人 = al ADV) a (11) 


WW(a, 5) 小 波 变 化 系数 随 a,b 变化 ,可 作出 以 2 为 横 坐 标 ,a 为 纵 坐 标的 二 维 小 波 系数 取 (ac,， 2 ) 等 值 线 

图 。 而 正 负 小 波 系数 的 突变 点 则 为 突变 点 。 在 时 间 域 上 对 关于 a 的 所 有 小 波 系 数 的 平方 进行 积分 ,得 到 小 波 
方差 , 见 公式 (12) ,反映 了 波动 能 量 在 时 间 尺 度 的 分 布 情况 ,从 而 确定 时 间 序 列 下 的 主 周期 。 

var(a) = | lw(a,b) [db (12 ) 


2.2.$S 双 累 积 曲线 法 
双 累 积 曲线 是 研究 两 个 参数 间 的 关系 一 致 性 及 变化 的 常用 方法 。 通 过 在 直角 坐标 系 中 绘制 一 个 变量 的 
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连续 累积 值 与 男 一 个 变量 的 连续 累积 值 的 关系 线 , 分 析 获 得 两 参数 间 关 系 的 一 种 常用 的 分 析 方 法 ”。 其 斜 
率 的 变化 可 以 作为 变量 间 关 系 发 生变 化 的 指示 ,拐点 出 现 的 时 间 则 为 变量 间 关 系 发 生 明显 变化 的 年 份 ”。 
该 方法 是 目前 在 水 文 气象 要 素 一 致 性 或 长 期 演变 趋势 分 析 中 最 简单 直观 和 广泛 应 用 的 分 析 方 法 ”。 当 只 
有 降水 的 变化 而 无 其 他 因素 影响 时 ,降雨 -径流 双 累 积 曲线 应 为 一 条 直线 ,而 当 有 人 类 活动 等 其 他 因素 影响 
时 ,曲线 将 发 生 偶 移 。 因 此 ,可 根据 曲线 发 生 偶 移 的 年 代 确 定 下 垫 面 受 人 类 活动 影响 的 时 间 ,而 根据 偶 移 程度 
可 推 央 人 类 活动 对 径流 影响 的 剧烈 程度 。 因 此 ,降雨 -径流 的 双 累 积 曲 线 可 以 反应 人 类 活动 对 径流 影响 的 阶 
段 性 特征 。 

双 累 积 曲 线 具 体 的 计算 过 程 如 下 : 

(1) 首 先 将 双 累 积 曲线 斜率 未 发 生 偏 移 的 时 期 作为 未 受 人 类 活动 干扰 的 基础 期 ,建立 基准 期 降雨 量 -2 了 


和 累积 径流 量 》 0 的 相关 方程 (13): 


>P=F》>O+DB (13 ) 
(2) 在 假设 变化 期 与 基准 期 的 气 修 变化 同步 的 前 提 下 ,公式 (14) 得 到 的 变化 期 模拟 径流 与 实测 径流 的 差 
值 即 为 径流 受 人 类 活动 影响 的 变化 量 
4AO" = mean( 0Q,) - mean( 0,) (14) 
(3) 变 化 期 实测 径流 与 基准 期 实测 径流 的 差 值 为 径流 受气 候 变 化 与 人 类 活动 共同 影响 的 变化 量 由 公式 
(15) 可 得 


AQ = mean( (0,,) - mhéantQ,,) (15) 
(14) (15) 式 中 ,下 标 cwm 分 别 代表 模拟 值 和 实测 值 ,1,2 分 别 代表 基准 期 和 变化 期 。 


3 结果 与 分 析 


3.1 径流 量 的 年 内 和 年 际 分 配 特征 
3.1.1 年 内 不 均匀 系数 

图 3 展示 了 密云 水 库 潮 河 与 白河 的 水 文 监测 站 点 1960 一 2014 年 的 人 库 平均 径流 流量 。 张 家 坟 和 下 会 水 
文 站 径流 量 的 季节 性 变化 特征 一 致 ; 主 汛期 7 一 9 月 份 径流 量 较 大 ,在 12.3 一 39.1m/s 之 间 ; 而 其 他 月 份 径流 
量 则 在 1.3 一 12.0m /s 之 间 。 两 个 站 作 比 对 ,可 见 张 家 坟 站 的 月 平均 流量 为 8.28m /s 是 下 会 站 的 月 平均 流量 
为 5.71m’/s 的 1.45 售 。 


利用 不 均匀 系数 可 以 有 效 反应 径流 量 的 年 内 分 配 es 
特征 。 图 4 展示 了 张 家 攻 站 和 下 会 站 不 均匀 系数 计算 80.00 


口 下 会 站 
国 张 家 坟 站 


结果 。 张 家 坟 站 不 均匀 这 数 波动 性 较 大 ,在 0.20 一 1.80 


E 
之 间 , 说 明 该 站 径流 量 的 年 内 变化 较 大 ,而 下 会 站 径流 三 昌 5000 
的 不 均匀 系数 在 0.04 一 0.68 之 间 。 潮 河流 域 的 人 库 径 “号 3000 
流年 内 变化 比 百 河流 域 人 库 径流 小 ,不 均匀 系数 显著 相 ”1000 
关 ( P<0.01 0 


去 -10.00 6 7 8 9 10 11 12 
3.1:2 年 际 变 化 -20.00 月 份 Month 
利用 5 年 滑动 平均 3 分别 对 下 会 站 和 张 家 坟 站 径 
三 一 、 上 一 a 家 坟 和 下 会 站 月 平均 流量 
流量 的 年 际 变 化 进行 分 析 , 如 图 5、 图 6 所 示 。50 年 来 ， Re ee dis 


Fig.3 The monthly runoff of Zhang Jiafen and Xiahui station 


白河 流域 和 潮 河 流域 人 库 量 在 线性 趋势 上 上 有 减少 的 

趋势 ,而 在 20 世纪 70 年 代 和 20 世纪 90 年 代 的 径流 量 有 增加 的 趋势 ,20 世纪 60 年 代 、80 年 代 和 2000 年 以 来 
径流 量 有 减少 趋势 。 其 中 从 图 5、 图 6 中 可 得 , 潮 河 流域 和 白河 流域 的 人 库 径 流量 在 1973 年 .1974 年 达到 最 
大 值 ,分 别 为 28.1.42.4m /s。 而 最 小 值 两 个 水 文 站 发 生 在 2002 年 ,其 中 下 会 站 为 1.86m /s, 张 家 坟 站 为 
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2.2m”/s。 密 云 水 库 50 年 来 入 库 径 流量 


1.8 


1.5 


不 均匀 系数 
Nonuniform coefficient 
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图 4 1960 一 2014 年 张 家 坟 和 下 会 站 不 均匀 系数 


Fig.4 Unevenness of annual runoff distribution of Zhang Jiafen and Xiahui Station from 1960 to 2014 
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图 5 潮 河 流域 年 径流 量变 化 及 $ 年 滑动 平均 


Fig.s “Runoff variation of Chao River basin and Sa Moving-Average 
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图 6 白河 流域 年 径流 量变 化 及 5 年 滑动 平均 


Fig.6 Runoff variation of Bai River basin and Sa Moving-Average 
32 扼要 分 性 


利用 Mann-Kendall 法 检验 分 析 了 在 显著 性 水 平 w = 0.05 下 密云 水 库 和 人 库 的 水 文 站 点 张 家 坟 站 (白河 ) 和 
下 会 站 ( 漳 河 )1960 一 2014 年 径流 变化 趋势 (图 7, 图 8)。 张 家 坟 站 和 下 会 站 径流 分 别 在 20 世纪 60 年 代 初 经 
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历 短暂 的 下 降 , 然 后 进入 径流 量 增加 的 阶段 ,该 过 程 一 直 持 续 到 20 世纪 80 年 代 , 随 后 径流 量 急剧 下 降 , 其 中 


张 家 坟 站 在 1987 年 发 生 突变 ,下 会 站 则 在 在 1998 年 发 生 突 变 。 


2 2 
UF 1 
8 0 
2 5 
泽 号 
有 “3 
曙 吊 
二 全 
RS NV 2 
-3 
-4 一 4 
1960 1970 1980 1990 “ 2000 2010 2020 1960 1970 1980 1990 “2000 | 2010 2020 
年 份 Year 年 份 -Year 
图 7 张 家 坟 站 径流 Mann-Kendall 突变 分 析 图 8 下 会 站 径流 Mann=:Kenidall 突变 分 析 


Fig.7 The Mann-Kendall abrupt change test on runoff of Zhang Fig.8 The Mann-Kendall abrupt change test on runoff of Xiahui 
Jiafen Station Station 
UF :表示 M-K 顺序 统计 值 ; UB: 表示 M-K 逆序 统计 值 UF :表示 M-K 顺序 统计 值 ; UB 表示 M-K 逆序 统计 值 


3.3 ”周期 性 变化 

采用 Morlet 小 波 对 密云 水 库 白河 和 漠河 流域 50 年 来 年 径流 量 进行 周期 演变 分 析 。 分 别 得 到 白河 张 家 坟 
站 漠河 下 会 站 小 波 变 换 系 数 实 部 图 (图 9, 图 11) 和 小 波 方 差 图 (图 10, 图 12)。 小 波 系数 的 变化 特征 可 以 用 
来 表征 年 径流 量 的 变化 特征 。 当 小 波 系 数 实 部 为 正 值 时 ,代表 丰 水 期 , 负 值 代表 枯水期 ,而 小 波 系数 实 部 为 0 
时 , 则 表示 年 径流 从 丰 水 期 转 为 枯水期 或 由 枯水期 转 为 丰 水 期 的 转折 点 。 

由 图 9 可 见 ,白河 流域 年 径流 量变 化 主要 存在 5 一 8a ,14 一 18a 周期 振荡 。 在 5 一 8a 的 特征 尺度 下 , 丰 枯 
交替 特别 频繁 ,并 以 7a 为 振荡 中 心 。 主 要 集中 在 20 世纪 60 年 代 一 20 世纪 90 年 代 之 间 ,其 中 20 世纪 60 年 
代 一 20 世纪 80 年 代 经 历 了 强烈 的 在- 枯 交 蔡 转 变 过 程 。14 一 18a 特征 尺度 下 ,主要 集中 在 20 世纪 90 年 代 以 
后 ,震荡 中 心 在 17a 左右 。 而 通过 图 10 ,白河 流域 存在 7a 和 14a 两 个 主 周期 。 其 中 在 7a 的 小 波 方差 值 最 大 ， 
是 第 一 主 周期 。 图 11 是 潮 河 流域 年 径流 量变 化 情况 ,主要 存在 $ 一 8a 的 周期 振荡 ,14 一 18a 周期 振荡 不 明显 。 
与 白河 流域 年 径流 量变 化 特征 类 和 似 , 在 5 一 8a 存在 丰 枯 频繁 交替 的 现象 ,以 7a 为 振荡 中 心 ,在 20 世纪 60 年 
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图 9 白河 流域 张 家 坟 站 小 波 变 换 系 数 实 部 图 


Fig.9 The contours of the real part of Morlet wavelet coefficients of runoff in Bai River basin 
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代 一 20 世纪 70 年 代 经 历 了 对 比 强烈 的 丰 - 村 交替 转变 过 程 。 而 结合 日 河流 域 年 径流 量 小 波 方 差 图 ( 图 10)， 
只 存在 7a 的 主 周期 ,14a 峰值 并 不 明显 。 因 此 认为 ,日 河流 域 年 径流 量 第 一 主 周期 为 7a。 因 此 ,对 于 白河 和 
漳 河 两 个 密云 水 库 主 要 入 库 河 流 来 说 ,径流 变化 趋势 赋 有 类 似 之 处 同时 也 存在 着 差异 。 


3.4 ”密云 水 库 入 库 径 流量 变化 驱动 因子 分 析 


0.7 

影响 径流 量 的 因素 主要 包含 两 个 方面 ,一 是 自然 因 i 

素 ,二 是 人 类 活动 。 自 然 因 素 主要 包含 地 形 、 土 壤 等 下 8 0 

垫 面条 件 和 气候 条 件 ,一般 认为 下 垫 面 在 一 定时 间 尺 度 共计 04 

上 是 不 发 生 改 变 的 ,因此 自然 因素 主要 考虑 的 是 气候 因 次 号 03 

素 如 降雨 对 径流 量 的 影响 。 人 类 活动 主要 包括 取 用 水 、 02 

水 利 工程 和 土地 利用 等 方面 。 因 此 ,采用 了 降水 - 径 0 
流 双 累积 曲线 对 降雨 和 人 类 活动 对 密云 水 库 人 库 径 流 0 5 10 15 l/ss VM)30 35 


量 的 影响 进行 分 析 。 

图 13 为 白河 流域 降水 -径流 双 累 积 曲 线 。 根 据 径 
流 -降雨 双 累 积 曲 线 的 斜率 发 生变 化 的 位 置 确定 1973 
年 和 1996 年 为 突变 点 ,将 白河 流域 径流 序列 分 为 
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图 10 自 河流 域 张 家 坟 站 小 波 方差 


Fig.10 The wavelet variances of annual runoff in Bai River basin 
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图 11 


潮 河 流域 下 会 站 小 波 变换 系数 实 部 图 


Fig.11 The contours of the real part of Morlet wavelet coefficients of annual runoff in Chao River basin 


1960 一 1973 年 , 1974 一 1996 年 , 1997==-2014 年 3 个 阶 


1.4 
段 。 其 中 ,第 一 阶段 41960 一 1973 年 可 以 作为 白河 流域 。12 
径流 序列 的 基准 期 -从 二 阶段 的 降水 -径流 双 黑 积 曲 线 染 & 
的 拟 合 方程 中 , 双 黑 积 的 县 线 斜率 从 基准 期 的 0.871 下 党 
降 到 0.3007 衰减 率 为 65.5% ,表明 人 类 活动 对 白河 流域 04 
的 径流 天 本 极 大 的 影响 。 图 14 是 潮 河流 域 降雨 -径流 0 
双 依 积 曲 线 ; | 同 百 河流 域 径 流 趋势 类 似 , 潮 河流 域 下 会 0 


站 的 径流 曲线 突变 点 为 1972 年 和 1997 年 ,将 第 一 阶段 
1960 一 1971 年 作为 基准 期 ,其 他 两 个 阶段 表示 了 人 类 
活动 对 径流 的 影响 。 通 过 对 双 累 积 曲 线 方 程 进行 线性 
拟 合 ,得 到 各 阶段 双 累 积 曲 线 的 斜率 , 潮 河 流域 基准 期 
的 斜率 为 0.4443 而 第 三 阶段 的 斜率 则 下 降 至 0.1469， 


时 间 尺 度 Time scales/a 


图 12 潮 河流 域 下 会 站 小 波 方差 
Fig.12 The wavelet variances of annual runoff in Chao 


了 River basin 


衰减 率 为 66.9% 。 可 见 人 类 活动 对 潮 河 流域 径流 同样 有 巨大 的 影响 。 表 1 列 出 了 人 类 活动 及 降雨 对 密云 水 
库 白 河和 潮 河 流域 径流 量 的 影响 。 表 明 潮 河流 域 和 白河 流域 径流 量 有 减少 趋势 。 昌 河流 域 人 库 流 量 从 
1960 一 1973 年 的 5.98 亿 mi ,减少 到 1997 一 2014 年 的 2.07 亿 mi ,减少 了 65.38%。 而 潮 河 流域 人 库 流量 则 从 


http://www.ecologica.cn 


6 期 秦 丽 欢 ， 等 :密云 水 库 上 游 径 流 变化 趋势 及 影响 因素 1949 


3.2 亿 m3(1960 一 1971 年 ) 减 少 到 1.19 亿 ms3(1998 一 2014 年 ) 减 少 了 62.81%。 白河 流域 人 库 径 流量 减少 幅度 
更 为 明显 。 通 过 对 潮 白 河流 域 突变 期 的 分 析 , 分 别 得 到 了 降雨 和 人 类 活动 对 入 库 人 径流 的 影响 比例 ( 表 1) 。 降 
雨 和 人 类 活动 对 白河 入 库 两 个 突变 期 的 径流 的 影响 在 8.21% 、10.10% 和 91.90% .89.90% ; 对 潮 河 流域 的 影响 
分 别 在 22.60% .27.64% 和 77.40% .72.36% 。 因 此 ,人 类 活动 可 能 是 导致 密云 水 库 人 库 径 流 减 少 的 主要 因素 。 
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图 13 ”白河 流域 张 家 坟 站 双 累 积 径流 曲线 图 14 潮 兆 认 域 下 朗 好 冯 累 稳 醋 流 曲 线 图 


Fig.13 “Double mass curve of the precipitation and runoff in Bai Fig.14 Double mass curve of the precipitation and runoff in Chao 


River River 
表 1 ， 人 类 活动 及 降水 变化 对 密云 水 库 入 库 径 流量 的 影响 ( 10%m”) 
Table 1 Influence of precipitation and human activities on runoff decline in the upper reaches of the Miyun River Basin 
水 文 乌 实测 英 拟 降雨 Precipitation 人 类 活动 Human activities 
流域 0 同 本 克 
人 ie easure iiio 变化 量 比例 /% 变化 量 比例 /% 
Stalon value value Change value Ratio Change value Ratio 
白河 张 家 攻 1960 一 1973 598.73 598.73 
1974- 一 1996 430.32 584.91 —13.82 8.21 —154.59 91.79 
1997 一 2014 207.35 546.79 -38.12 10.10 —339.44 89.90 
潮 河 下 会 1960 一 1971 320.72 320.72 
1972 一 1997 209.17 295.51 =25.21 22.60% —86.34 77.40 
1998 一 2014 118.93 264.95 -55.77 27.64% —146.02 72.36 
4 讨论 


4.1 降雨 对 径流 影响 分 析 

密云 水 库 白 河流 域 和 少 河 流域 径流 具有 年 内 分 配 不 均匀 的 特性 ,在 年 际 变化 中 径流 量 有 缓慢 减 少 的 趋 
势 。 高 迎春 等 “研究 以 长 时 间 序 列 降 十 和 径流 的 数据 进行 分 析 , 发 现 随 着 时 间 的 向 后 推移 ,降雨 与 径流 之 间 
的 关系 逐渐 减弱 。 而 降雨 在 相同 时 间 尺 度 上 并 没有 发 生 明 显 的 变化 ,降雨 对 径流 的 影响 比例 变化 不 大 (如 白 
河流 域 大 第 二 阶段 的 8.21% 变 化 为 第 三 阶段 的 10.10% ) ( 表 1) ,从 这 一 方面 可 以 推断 ,降雨 可 能 不 是 影响 径 
流量 变化 的 关键 因素 。 根 据 表 1 中 结果 ,降雨 对 潮 河 流域 径流 减少 的 影响 在 1974 一 1996 年 和 1997 一 2014 年 
两 个 阶段 分 别 为 22.60% 和 27.64% 。 李 子 君 和 李 秀 彬 ”通过 分 析 降 十 和 径流 变化 的 一 致 性 也 得 到 密云 水 库 
潮 河流 域 径 流 减少 受降 雨 影 响 较 小 的 结论 。 而 对 降雨 对 白河 流域 的 影响 1974-1996 年 和 1997 一 2014 年 两 个 
阶段 分 别 为 8.21% 和 10.10%。 王 静 怡 等 .同样 通过 分 析 入 库 径 流 和 降雨 的 关系 中 进一步 断定 ,降雨 不 是 决 
定 径 流 变 化 的 单一 因素 。 王 巧 平等 ”在 海河 流域 研究 了 人 类 活动 对 径流 的 影响 ,结果 表明 ,降雨 能 够 影响 该 
流域 径流 量 的 减少 ,但 不 是 主要 因素 。 对 海河 流域 而 言 ,降雨 主要 集中 在 7、8 月 份 , 产 流 也 主要 由 汛期 的 几 次 
大 的 降雨 形成 ,影响 产 流 的 因素 主要 来 自 于 土壤 吸水 .植物 蒸 散 发 , 填 洼 等 损失 。 由 于 雨 期 的 蒸发 量 可 以 忽略 
不 计 , 而 下 热 面 的 变化 对 降雨 的 影响 可 以 忽略 不 计 , 因 此 ,地 表 径 流 的 变化 主要 来 自 人 类 活动 引起 地 表 植 被 改 
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变 蒸腾 作用 的 等 方面 的 影响 。 
4.2 ”人 类 活动 对 径流 影响 分 析 

人 类 活动 如 土地 利用 变化 ,水 库 的 拦截 以 及 器 流域 的 水 库 补 水 是 导致 径流 变化 的 主要 原因 。 如 表 1 中 所 
示 , 人 类 活动 对 白河 流域 潮 河 流域 影响 分 别 在 89.90% 一 91.79% 和 72.36% 一 77.40%。 通 过 对 流域 内 赤峰 县 、 
丰 宁 县 和 洪 平 县 的 人 类 活动 用 水 量 进行 统计 ,发 现 1985 一 2015 年 间 用 水 量 呈 增加 趋势 ,总 用 水 量 从 0.02 亿 
mm 增加 到 0.14 亿 m 一 。 因 此 流域 用 水 量 的 增加 ,导致 潮 白 河 径流 量 的 减少 。 同 时 ,密云 水 库 流 域 从 1990 一 
2010 年 间 ,土地 利用 类 型 发 生 了 重大 的 转变 ,农田 .草地 和 水 体 的 面积 分 别 减 少 了 30% .48% .61% ,林地 增加 
了 30%。 由 于 林地 相 比 农田 和 草地 具有 较 强 的 水 分 蒸腾 作用 ,导致 径流 量 减少 。 而 小 流域 的 综合 治理 提 
高 了 植被 覆盖 率 ,增加 了 植被 的 截流 量 , 因 此 植物 截留 量 .土壤 持 水 量 及 蒸 散 发 量 等 的 增 大 ,使 得 径流 量 减 
少 ”” ,与 此 同时 ,密云 水 库 上 游 有 26 座 水 库 ,总 库容 达 1.325 亿 m*… ,水 库 会 在 一 定 程度 上 增加 东 发 损失 和 
渗 漏 量 ,从 而 减少 径流 量 。 除 此 之 外 ,流域 内 水 库 向 流域 外 调 水 ,也 在 一 定 程度 上 减少 径流 y 如 白河 流域 的 白 
河 堡 水 库 承 担 防 洪 灌溉 的 功能 外 ,还 向 流域 外 的 官厅 水 库 补 水 ,平均 每 年 补水 量 达 -0.78 亿 m55 回 十 三 陵 水 库 
平均 每 年 补水 0.1m”” ,这 样 向 流域 外 水 库 补 水 的 过 程 ,在 一 定 程度 上 减少 潮 白 河流 域 的 径流 量 。 
4.3 ”南水北调 重大 调 水 工程 影响 下 ,密云 水 库 水 量 调度 的 启示 

密云 水 库 在 南水北调 工程 中 起 到 了 对 进 京 水 量 调 蓄 的 作用 ,并 于 2015 年 9 月 开始 接纳 南 水 进 入 , 蓄 水 量 
得 到 了 明显 的 补给 ,水 位 一 度 上 升 至 140m。 南 水 北 调 之 水 进入 密 去 水库; 一 方面 缓解 了 水 库 的 供水 压力 ,有 
效 补 充 了 水 库 的 水 量 , 另 一 方面 也 为 水 库 运 行 带 来 风险 。 通 过 对 卡 游 径流 变化 特征 的 研究 ,将 其 与 南水北调 
入 库 调 蓄 流 量 结合 起 来 , 丰 水 年 的 汛期 结合 小 调 水 量 , 枯 水 年 非 讯 期 结合 大 调 水 量 。 这 样 可 以 有 效 协调 和 发 
挥 上 游 入 库 径 流 和 南水北调 之 水 的 特性 和 优势 ,保证 密云 水 库 正常 的 蓄 水 供水 功能 。 


5 结论 


潮 河 和 白河 作为 密云 水 库 上 游 两 条 主要 的 入 库 河流 ,对 密云 水 库 有 重要 的 影响 和 作用 。 本 文 分 别 对 潮 河 
流域 和 日 河流 域 径流 进行 了 分 析 比 较 ,结论 如 下 : 

(1) 从 50 年 径流 数据 中 可 以 看 出 , 济 河 与 日 河 的 径流 量 年 内 分 配 不 均匀 的 情况 ,主要 的 径流 量 多 集中 在 
汛期 (7 一 9 月 ) ,汛期 径流 量 占 全 年 比例 70%。 与 潮 河 流域 径流 量 相 比 ,白河 流域 的 年 内 变化 更 为 剧烈 ,径流 
量 减少 趋势 更 为 明显 。 

(2) 白 河 与 潮 河 流域 径流 部 具有 周期 性 的 变化 ,其 中 日 河流 域 径 流 存在 7a 和 17a 的 周期 ,而 潮 河 流域 径 
流 周期 为 7a。 

(3) 通 过 对 降 南 和 人 类 活动 对 密云 水 库 径 流量 的 影响 分 析 , 结 果 表 明白 河流 域 径流 受到 人 类 活动 变化 更 
为 剧烈 。 
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